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Abstract 

With its unique data samples at energies of 3.8-4.6 GeV, the BESIII experiment made a signif¬ 
icant contribution to the study of charmonium and charmonium-like states, i.e., the XYZ states. 
We review the results for observations of the Zc(3900) and Zc(4020) states, the A(3872) in e'''e“ 
annihilation, and charmonium i/'(l ^ 02 ) state, as well as measurements of the cross-sections of (jJXcJ 
and gJ/ij, and the search for e'’'e“ —)■ yXcj and 7y(4140). We also present data from BESIII that 
may further strengthen the study of the XYZ and conventional charmonium states, and discuss 
perspectives on future experiments. 
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I. INTRODUCTION 


Many charmonium and charnioniuni-like states were discovered at 5-factories in the first 
decade of the 21st century [l|. Whereas some of these are good charmonium candidates, 
as predicted in different models, many states have exotic properties, which may indicate 
that exotic states, such as multi-quark, molecule, hybrid, or hadron-quarkonium, have been 
observed [^. 

The BESIII experiment at the BEPCII storage ring started its first collisions in the 
tau-charm energy region in 2008. After a few years running at suitable energies for its well- 
defined physics programs j^, i.e., at J/'0 and peaks in 2009 and the '0(3770) peak 

in 2010 and 2011, the BESIII experiment started to collect data for the study of the XYZ 
particles, which were not described in the Yellow Book j^. 

As the design center-of-mass (c.m.) energy of the BEPCII was 2.0-4.2 GeV, there were 
not many options for data samples relevant to the XYZ-related physics. BESIII took its 
hrst data sample at the peak of "0(4040) in May 2011, with the aim of searching for the 
well-known Y(3872) in the 0(4040) radiative transition and possibly the excited P-wave 
charmonium spin-triplet states in similar transitions. This sample is about 0.5 fb“^, which 
is limited by the one-month running time left after the 0(3770) data taken in the 2010-2011 
run. Data were not collected at the 0(4170) peak because the CLEO-c experiment had 
already collected a sample of about 0.6 fb“^ for the study of Dg decays, which could be used 
for similar studies. 

The upgrade of BEPCIPs LINAC in summer 2012 increased the highest beam energy 
from 2.1 to 2.3 GeV, making it possible to collect data at higher c.m. energies (up to 
4.6 GeV). This made data collection possible at almost all known vector states, including 
Y(4260), Y(4360), 0(4415), and (marginally) Y(4660). 

The data collected at a c.m. energy of ^/s = 4.26 GeV turned out to be very fruitful. 
One month’s data of 525 pb“^ (from December 14, 2012 to January 14, 2013) produced 
observations of the charged charmonium-like state Zc(3900) j^, resulting in changes to the 
data collection plan for the 2012-2013 run. More data were accumulated at c.m. energies 
of 4.26 GeV and then 4.23 GeV. Data from the Y(4360) peak were also obtained in spring 
2013, and data from even higher energies (4.42 and 4.6 GeV) were recorded in 2014 after 
a hne scan of the total hadronic cross-sections between 3.8 and 4.6 GeV at more than 100 
energy points, with a total integrated luminosity of about 800 pb“0 

The dedicated data samples for the XYZ study are presented in Table [H which lists the 
nominal c.m. energy, measured c.m. energy, and integrated luminosity at each energy point. 
These data were used for all the analyses presented in this article. 


II. CHARGED CHARMONIUM-LIKE STATES: Z^S 

The BESIII experiment observed, for the first time, a charged charmonium-like state close 
to the DD* threshold Zc(3900)/Zc(3885), and a charged charmonium-like state close to the 
D*D* threshold Zc(4020)/Zc(4025) . Their neutral partners were also observed, confirming 
their isospins to be one. 
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TABLE I: The measured c.m. energy integrated luminosity [7| of each data sample collected 
for the study of the XYZ states. The uncertainties on the integrated luminosities are statistical 
only; a 1% systematic uncertainty common to all the data points is not listed. 


Data sample 

c.m. energy (MeV) 

C (pb-i) 

3810 

3807.65±0.10±0.58 

50.54±0.03 

3900 

3896.24±0.11±0.72 

52.61±0.03 

4009 

4007.62±0.05±0.66 

481.96±0.01 

4090 

4085.45±0.14±0.66 

52.63±0.03 

4190 

4188.59±0.15±0.68 

43.09±0.03 

4210 

4207.73±0.14±0.61 

54.55±0.03 

4220 

4217.13±0.14±0.67 

54.13±0.03 

4230 

4226.26±0.04±0.65 

1091.74±0.15 

4245 

4241.66±0.12±0.73 

55.59±0.04 

4260 

4257.97±0.04±0.66 

825.67±0.13 

4310 

4307.89±0.17±0.63 

44.90±0.03 

4360 

4358.26±0.05±0.62 

539.84±0.10 

4390 

4387.40±0.17±0.65 

55.18±0.04 

4420 

4415.58±0.04±0.72 

1073.56±0.14 

4470 

4467.06±0.11±0.73 

109.94±0.04 

4530 

4527.14±0.11±0.72 

109.98±0.04 

4575 

4574.50±0.18±0.70 

47.67±0.03 

4600 

4599.53±0.07±0.74 

566.93±0.11 


A. Observation of Zc(3900) and Zc(3885) 

1. Observation of Zc{3900) 


The BESIII experiment studied the e’''e“ —?■ 7r~^7r~ J/fj process at a c.m. energy of 
4.26 GeV using a 525 pb“^ data sample [^. A structure at around 3.9 GeV/c^ was observed 
in the TT^J/rf mass spectrum with a statistical significance larger than Scr, which is referred 
to as the Zc(3900). A fit to the 7r^J/'0 invariant mass spectrum (see Fig. [T]), neglecting 
interference, results in a mass of (3899.0±3.6±4.9) MeV/c^ and a width of (46±10±20) MeV. 
The associated production ratio was measured to be i? = = 

(21.5 ±3.3 ±7.5)%. 


In the Belle experiments, the cross-section of —?■ Jjif was measured from 

3.8-5.5 GeV using the initial state radiation (ISR) method. The intermediate states in 
V(4260) —decays were also investigated [^. The Zc(3900) state (referred to as 
Z(3900)’'" in the Belle paper) with a mass of (3894.5 ± 6.6 ± 4.5) MeV/c^ and a width of 
(63 ±24 ±26) MeV was observed in the n^J/ip mass spectrum (see Fig. [T]) with a statistical 
significance larger than 5.2cr. 


The %c(3900) state was confirmed shortly after with GLFO-c data at a c.m. energy of 
4.17 GeV [91], and the mass and width agreed very well with the BFSIII and Belle measure¬ 
ments. 


A neutral state Zc(3900)‘^ —)■ with a significance of 10.4cr was observed at BESIII 
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FIG. 1: Unbinned maximum likelihood fit to the distribution of the M^i^^TrJ/'ijj) (left panel from 
BESIII and right panel from Belle). Points with error bars are data, the curves are the best 
ht, the dashed histograms are the phase space distributions and the shaded histograms are the 
non-7r+7r“ J/'i/’ background estimated from the normalized J/V' sidebands. 


in e^e“ —)■ /ip with c.m. energy ranges from 4.19-4.42 GeV jJ. The mass and width 

were measured to be (3894.8 ± 2.3 ± 3.2) MeV/c^ and (29.6 ± 8.2 ± 8.2) MeV, respectively. 
This state is interpreted as the neutral partner of the Zc(3900)^, as it decays to /ip and 
its mass is close to that of Zc(3900)^. This is in agreement with the previously reported 
3.5cr evidence for Zc(3900)° in the CLEO-c data i. The measured Born cross-sections of 
e’''e“ —)■ 7r°7r°J/'0 were about half of those for —)■ 'K~^n~J/ip measured in the Belle 

experiment j^, which is consistent with the isospin symmetry expectation. 


2. Observation of Zc{3885) 


The Zc(3900) state observed in the nJ/ip hnal state is close to and above the DD* mass 
threshold. With the same data sample at ^/s = 4.26 GeV, the BESIII experiment studied 
e'''e“ —)■ 7r^(ZlZ)*)^. A structure (referred to as Zc(3885)) was observed in the [DD*)^ 
invariant mass distribution When htted to a mass-dependent-width Breit-Wigner 

(BW) function, the pole mass and width were determined to be (3883.9± 1.5±4.2) MeV/c^ 
and (24.8 ± 3.3 ± 11.0) MeV, respectively (see Fig. [S]). The angular distribution of the 
Zc(3885) system favors a = I’*' assignment for the structure and disfavors 1“ or 0“. The 
production rate was measured to be cr(e'''e“ 7r^Zc(3885)^) x i3(Zc(3885)^ {DD* 

(83.5 ±6.6 ±22.0) pb. 




An important question is whether Zc(3885) is the same as Zc(3900) ji, Isf. The mass 
and width of Zc(3885) are 2a and la, respectively, below those of Zc(3900), as observed 
by the BESIII and Belle experiments. However, neither £t considers the possibility of 
interference with a coherent non-resonant background, which could shift the results. A spin- 
parity quantum number determination for Zc(3900) would provide an additional test of this 
possibility. 


Assuming the Zc(3885) structure is caused by Zc(3900), we obtain pl^z‘'( 39 ooj^!^j/J) ~ 
6.2 ± 1.1 ± 2.7. This ratio is much smaller than typical values for decays of conventional 
charmonium states above the open charm threshold, e.g., r('0(377O) —> DD)/r(ip(3770) 
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FIG. 2: The M{D^D*-) (left) and M{D+D*^) 
4.26 GeV. The curves show the best fits. 
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(right) distributions for selected events at y/s = 


TT+TT-J/ij) = 482 ±84 [l2| and r(V'(4040) ^ ±)W±)W)/r(V'(4040) ^ r]J/ij) = 192 ±27 |l3 


This suggests very different dynamics in the Y (4260)-Zc(3900) system. 


B. Observation of .^c(4020) and .Zc(4025) 


1. Observation o/Zc(4020) 


BESIII measured e'''e“ TT^TT~hc cross-sections [1^ at c.m. energies of 3.90-4.42 GeV. 
Intermediate states were studied by examining the Dalitz plot of the selected 7r~^7r~hc can¬ 
didate events. The he signal was selected using 3.518 < M^rje < 3.538 GeV/c^, and 7r+7r“hc 
samples of 859 events at 4.23 GeV, 586 events at 4.26 GeV, and 469 events at 4.36 GeV were 
obtained with purities of 65%. Although there are no clear structures in the 7r’''7r“ system, 
there is clear evidence for an exotic charmonium-like structure in the n^hc system, as clearly 
shown in the Dalitz plot. Figure |3] (left) shows the projection of the M{n^hc) (two entries 
per event) distribution for the signal events, as well as the background events estimated 
from normalized mass sidebands. There is a signiheant peak at around 4.02 GeV/c^ 
(Zc(4020)), and there are also some events at around 3.9 GeV/c^ (inset of Fig. E] (left)), 
which could be Zc(3900). The individual datasets at y/s = 4.23, 4.26, and 4.36 GeV show 
similar structures. 


An unbinned maximum likelihood £t was applied to the M{7r^hc) distribution summed 
over the 16 rje decay modes. The data at 4.23, 4.26, and 4.36 GeV were htted simultaneously 
to the same signal function with common mass and width. Figure E] (left) shows the htted 
results. The mass and width of Zc(4020) were measured to be (4022.9 ± 0.8 ± 2.7) MeV/c^ 
and (7.9 ± 2.7 ± 2.6) MeV, respectively. The statistical signihcance of the Zc(4020) signal 
was found to be greater than 8.9cr. 

Adding Zc(3900) with mass and width hxed to the BFSIII measurements to the 
£t results in a statistical signihcance of 2.la (see the inset of Fig. El (left)). At the 90% 
conhdence level (G.L.), the upper limits on the production cross-sections are set to a{e'^e~ 
7r^Zc(3900)^ —)■ 7i^7i~hc) < 13 pb at 4.23 GeV and < 11 pb at 4.26 GeV. These are lower 
than those of Zc(3900) j^. 
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FIG. 3: Sum of the simultaneous fits to the M{'K^hc) (left panel) and M{'K^hc) (right panel) 
distributions at 4.23, 4.26, and 4.36 GeV in the BESIII data; the inset in the left panel shows the 
sum of the simultaneous fit to the distributions at 4.23 and 4.26 GeV with Zc(3900) and 

Zc(4020). Dots with error bars are data; shaded histograms are normalized sideband background; 
the solid curves show the total fit, and the dotted curves the backgrounds from the fit. 


e'''e TT^Ti^hc at y/s = 4.23, 4.26, and 4.36 GeV for the first 
time ll5i. The measured Born cross-sections were about half of those for e’''e“ Tr+vr^hr, 


BESIII also observed 


which agree with expectations based on isospin symmetry within systematic uncertainties. 
A narrow structure with a mass of (4023.9 ± 2.2 ± 3.8) MeV/c^ (for htting, the width was 


hxed to that measured in the e+e“ 7r^7r~hc process [l^ because of low statistics) was 
observed in the Ti^hc mass spectrum (Fig. E] (right)). This structure is most likely the neutral 
isospin partner of the charged Zc(4020) observed in the e^e“ —)■ 7i~^n~hc process [l^. This 
observation indicates that there are no anomalously large isospin violations in vTTrhc and 
7 rZc( 4020 ) systems. 


2. Observation of Zc{‘i025) 


The BESIII experiment also studied the e+e —?■ {D*D*)^7r^ process at 4.26 GeV using 


an 827 pb“^ data sample [16[. Based on a partial reconstruction technique, the Born cross- 
section was measured to be (137 ± 9 ± 15) pb. A structure near the {D*D*)^ threshold in 
the TT^ recoil mass spectrum was observed, and this is denoted as Zc(4025) (see Fig. 0] (left)). 
The measured mass and width of the structure were (4026.3 ±2.6 ±3.7) MeV/c^ and (24.8 ± 
5.6 ± 7.7) MeV, respectively, from the £t with a constant-width BW function for the signal. 
The associated production ratio 
0.09 ±0.06. 


a{e+e (4025K^(g*5*) determined to be 0.65 ± 

a(e^e —¥(D*D*)^7r+) 


Using data at ^/s = 4.23 and 4.26 GeV, a structure was observed in the recoil mass 
spectrum in the e’''e“ —)■ D*^D*^{D*~^D*~)7i^ process 1^. Assuming that the enhancement 
is due to a neutral state decaying to D*D*, the mass and width of its pole position were 
determined to be (4025.5t|;7±3.1) MeV/c^ and T = (23.0±6.0±1.0) MeV, respectively (see 
Fig. m (right)). The Born cross-section a{e'^e~ —)■ Zc(4025)°7r° (^D*ojj*o _|_ /}*+£)*-) 7 j- 0 ) 
was measured to be (61.6 ± 8.2 ± 9.0) pb at 4.23 GeV and (43.4 ± 8.0 ± 5.4) pb at 4.26 GeV. 
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FIG. 4: Unbinned maximum likelihood fit to the vr^ recoil mass spectrum (left) in e'’'e“ —^ 
[D *at ^/s = 4.26 GeV, and to the vr^ recoil mass spectrum (right) in e"'“e“ —^ {D*D*)^tt^ 
at y/s = 4.23 and 4.26 GeV at BESIII. 


The ratio ^^( 4025 U J-^(d*d*) +^-) compatible with unity at y/s = 4.26 GeV, which is 

expected from isospin symmetry. In addition, Zc(4025)° has a mass and width that are very 
close to those of Zc(4025)^, which couples to {D*D*)^. Therefore, the observed Zc(4025)° 
state is a good candidate for the isospin partner of Zc(4025)^. 

As the Zc(4025) parameters agree to within 1.5cr with those of Zc(4020), it is very probable 
that they are the same state. As the results for Zc(4025)^ are only from data at 4.26 GeV, 
extending the analysis to 4.23 GeV and 4.36 GeV will probably provide a dehnite answer. 


III. OBSERVATION OF V(4260) ^ 7X(3872) 

BESIII observed 7A(3872) 77r+7r~JM, with J/'ip reconstructed through its 

decays into lepton pairs = e^e~ or /i^/i“) [l8|. 

The M'K~ J/Ip) distribution (summed over all energy points), as shown in Fig. [5] (left), 
was htted to extract the mass and signal yield of X(3872). The ISR 'ip{2S) signal was used 
to calibrate the absolute mass scale and to extract the resolution difference between the data 
and a Monte Garlo (MG) simulation. Figure |5] shows the htted result: the measured mass 
of X(3872) was (3871.9 ± 0.7 ± 0.2) MeV/c^. From a ht with a hoating width, we obtain a 
width of (O.Olg o) MeV, or less than 2.4 MeV at the 90% G.L. The statistical signihcance of 
A:(3872) is 6.3a. 

The Born-order cross-section was measured, and the results are listed in Table m For 
4.009 and 4.36 GeV data, because the X(3872) signal is not signihcant, upper limits on the 
production rates are given at the 90% G.L. The measured cross-sections at around 4.26 GeV 
are an order of magnitude higher than the NRQGD calculation of continuum production jl^ . 
which may suggest the X(3872) events come from resonance decays. 

The energy-dependent cross-sections were htted with a Y (4260) resonance (parameters 
hxed to PDG [1^ values), linear continuum, or Fl-transition phase space (oc F^) term. 
Figure |5] (right) shows all the htted results, which imply that y^/ndf = 0.49/3 (G.L. = 
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FIG. 5: Left panel: fit the M(tt~^ 7r~ J/ip) distribution observed at BESIII. Dots with error bars 
are data, the curves are the best fit. Right panel: fit to a^[e~^e~ —)■ 7X(3872)] x B[X(3872) — >• 
7:'^7r~ J/ip] measured by BESIII with a y(4260) resonance (red solid curve), a linear continuum 
(blue dashed curve), or an El-transition phase space term (red dotted-dashed curve). Dots with 
error bars are data. 


TABLE II: The product of the Born cross section (T'®(e+e —)• 7^(3872)) and B{X{3872) —)• 

TT~^Tr~ Jjijj) at different energy points. The upper limits are given at 90% C.L. 


y/s (GeV) a^[e^e 

^ 7A(3872)] • R(A(3872) ^ vr+vr-J/V>) (pb) 

4.009 

0.00 ±0.04 ±0.01 or < 0.11 

4.229 

0.27 ±0.09 ±0.02 

4.260 

0.33 ±0.12 ±0.02 

4.360 

0.11 ±0.09 ±0.01 or < 0.36 


92%), 5.5/2 (C.L. = 6%), and 8.7/3 (C.L. = 3%) for a F(4260) resonance, linear continuum, 
and phase space distribution, respectively. Thus, the Y (4260) resonance describes the data 
better than the other two options. 


These observations strongly support the existence of the radiative transition process 
y(4260) 7X(3872). The y(4260) 7^(3872) process could be another previously 

unseen decay mode of the Y (4260) resonance. Together with the transitions to the charged 
charmonium-like state Zc(3900) j5|,|8|,[^, this suggests that there might be some commonality 
in the nature of JA(3872), y(4260), and Zc(3900), and so the model developed to interpret 
any one of them should also consider the other two. As an example, the authors of Ref. 
integrated these states into a molecular picture to calculate e^e“ —)■ 7^(3872) cross-sections. 


Combining the above with the e'*'e —)■ tt+tt J/ip cross-section measurement at 

s = 4.26 GeV from BESIII 0, we obtain ^ 7A:(3872)] ■ B[X{m2) 

J/'0]/o"®(e+e“ —)■ TT+TT” J/-^) = (5.2 ±1.9) x 10“^, under the assumption that X(3872) 
and 7r+7r“J/'0 are only produced from y (4260) decays. If we take B[X(3872) —)■ 7r+7r“J/-^] = 
— then7^ = ^™^ 


7r'''7r 


5% 21 
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B(y(4260)^7r+7r-J/f/') 
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FIG. 6; Simultaneous fit to the Mi.ecoii(7r''"7r“) distribution of 7Xci events (left) and 7 Xc 2 events 
(right), respectively. Dots with error bars are data, red solid curves are total fit, dashed blue curves 
are background, and the green shaded histograms are J/ij: mass sideband events. 


IV. OBSERVATION OF ^(l^Da) 


BESIII observed A(3823) in the —)■ 7r’''7r“A(3823) —)■ 7r’''7r“7Xci process with a 

statistical signihcance of 6.2cr in data samples at c.m. energies of ^/s =4.23, 4.26, 4.36, 4.42, 
and 4.60 GeV j^. 

Figure [6] shows the htted results to tt+tt” recoil mass distributions for events in the Xci 
and Xc 2 signal regions. The ht yields 19 ± 5 A(3823) signal events in the yxci mode, with a 
measured mass of X(3823) of (3821.7 ± 1.3 ± 0.7) MeV/c^, where the hrst error is statistical 
and the second systematic. For the 7 Xc 2 mode, no signihcant A(3823) signal was observed, 
and an upper limit on its production rate could be determined. The limited statistics do not 
allow a measurement of the intrinsic width of A(3823). From a ht using the BW function 
(with a hoating width parameter) convolved with Gaussian resolution, it can be determined 
that r[A(3823)] < 16 MeV at the 90% G.L. (including systematic errors). This measurement 
agrees well with the values found by Belle [231]. 

The production cross-sections of cr^(e+e“ —)■ 7r+7r“A(3823)) ■ i3(A(3823) —)■ 7Xci, lXc 2 ) 
were also measured at these c.m. energies. The cross-sections of e’''e“ —)■ 7r’''7r“A(3823) were 
htted with the Y (4360) shape or the -0(4415) shape, with their resonance parameters hxed to 
the PDG values [1^. Figure 0 shows the htted results, which give the Kolmogorov-Smirnov 
statistic = 0.151 for the F(4360) hypothesis (HI) and = 0.169 for the -0(4415) 
hypothesis (H2), based on the Kolmogorov-Smirnov test [1^. Thus, both the K(4360) and 
-0(4415) hypotheses -^^obs < -^s.o.i = 0.509) are accepted at a 90% G.L. 

The A(3823) resonance is a good candidate for the 0(1 ^D 2 ) charmonium state. Accord¬ 
ing to potential models jl^, the H-wave charmonium states are expected to be within a mass 
range of 3.82-3.85 GeV. Among these, the 1 ^D 2 —)■ 7Xci transition is forbidden because of 
G-parity conservation, and the amplitude for 1 —)■ 'jXci is expected to be small jl^. The 

mass of 0(1 ^D 2 ) is in the 3.810 ~ 3.840 GeV/c^ range predicted by several phenomenolog¬ 
ical calculations [1^. In this case, the mass of 0(1 ^Zl 2 ) was above the DD threshold but 
below the DD* threshold. Because 0(1 ^Zl 2 ) DD violates parity, 0(1 ^Zl 2 ) is expected 
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FIG. 7: Comparison of the energy-dependent cross sections of cr^[e'''e“ —)• 7r+7r“X(3823)] • 
i3(X(3823) — IXci) to the y(4360) and ^(4415) line shapes. Dots with error bars (statistical 
only) are data. The red solid (blue dashed) curve shows a fit with the y(4360) (V’(4415)) line 
shape. 


to be narrow, in agreement with the observation, and '0(1 ^D 2 ) 7Xci is expected to be 
a dominant decay mode j^ . . From the cross-section measnrement, we obtain the ratio 

g[x( 3823 )^ 7 x‘'i] ^ (where systematic nncertainty cancels) at the 90% C.L., which also 
agrees with expectations for the '0(1 ^ 02 ) state 


28 


V. SEARCH FOR e+e” ^ 7Xcj AND e+e” -)■ 7^(4140) 


Using data samples collected at ^/s = 4.009, 4.23, 4.26, and 4.36 GeV, BESIII searched 
for IXcj {J = 0,1,2) with the snbseqnent decay XcJ 7 J /0 and J/0 

Evidence for the IXci and e'''e“ —)■ xXc2 processes was observed with a statistical 

signihcance of d.Ocr and 3Aa, respectively jl^. No evidence for XXco was observed. 

The corresponding Born cross-sections of —?■ xXcJ at different c.m. energies were 

calculated. Under the assumption that Xcj signals were absent, the upper limits on the 
Born cross-sections were calculated at the 90% C.L. These upper limits on the Born cross- 
section of e^e“ —?• xXcj are compatible with the theoretical prediction from the NRQCD 


calculation 19 


BESIII searched for Y (4140) via e’''e“ ^ 70<//'0 at -/s = 4.23, 4.26, and 4.36 GeV, but no 
significant y(4140) signal was observed in any of the data samples |30|. The upper limits of 
the product of the cross-section and branching fraction a[e'^e~ —?■ yV(4140)] ■ B{Y (4140) —?■ 
(pJ/tp) at the 90% C.L. were estimated to be 0.35, 0.28, and 0.33 pb at i/s = 4.23, 4.26, and 
4.36 GeV, respectively. 


These upper limits can be compared with the X(3872) production rates 18 


which 
-)■ 


were measured from the same data samples (see Sec. HID. The latter were ale~^e 
7^(3872)] ■ B(X(3872) ^ tt+tt-J/^) = [0.27 ± 0.09(stat) ± 0.02(syst)] pb, [0.33 ± 
0.12(stat) ± 0.02(syst)] pb, and [0.11 ± 0.09(stat) ± O.Ol(syst)] pb at = 4.23,4.26, and 
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4.36 GeV, respectively, which are of the same order of magnitude as the upper limits of 
(T[e+e“ —)■ 7^(4140)] ■ i3(F(4140) —)■ (pJ/ip) at the same energy. 

The branching fraction ;B(y(4140) —)■ (jyJ/'ilj) has not previously been measured. Using 
the partial width of y(4140) calculated under the molecule hypothesis j3l[ and the 

total width of the Y (4140) measured by CDF js^, the branching fraction was estimated to be 
approximately 30%. A rough estimation for B{X{3872) —)■ 7r^7r~ J/ip) is 5% jlH. Combining 
these numbers, the ratio a[e'^e~ 7^(4140)] /a—)■ 7X(3872)] can be estimated to be 
of the order of 0.1 or less at -y/s = 4.23 and 4.26 GeV. 


VI. STRUCTURES IN e+e —>■ charmonium + hadrons 


A. Observation of e+e —>■ uXcO 


Based on data samples collected between y/s = 4.21 and 4.42 GeV, the —)■ ooXco 

process was observed at y/s = 4.23 and 4.26 GeV for the hrst time (ssf, and the Born cross- 
sections were measured to be (55.4 ± 6.0 ± 5.9) and (23.7 ± 5.3 ± 3.5) pb, respectively. For 
other energy points, no signihcant signals were found, and upper limits on the cross-section 
at the 90% C.L. were determined. 


The data reveal a sizable uXco production at around 4.23 GeV/c^, as predicted in Ref. [34 
By assuming the uxco signals come from a single resonance, the TeeB{uXco), mass, and width 
of the resonance were determined to be (2.7 ± 0.5 ± 0.4) eV, (4230 ± 8 ± 6) MeV/c^, and 
(38 ± 12 ± 2) MeV, respectively (the ht is shown in Fig. [H]). The mrameters are consistent 
with those of the narrow structure in the Ti^7i~hc process [35|, but are not consistent 

with those obtained by htting a single resonance to the J/iIj cross-section (with 
a mass of (4259 ± 8lg) MeV/c^ and a width of (88 ± 23 ^ 4 ) MeV). This suggests that the 
observed ooxco signals are unlikely to originate from Y (4260). 


The uJXci ,2 channels were also sought, but no signihcant signals were observed; 

upper limits at the 90% C.L. on the production cross-sections were determined. The very 
small measured ratios of —)■ cross-sections to those for e+e” —)■ uxco me incon¬ 

sistent with the prediction in Ref. |37|. 


B. Measurement of e+e 


Usii^data samples collected at energies of 3.81-4.60 GeV, BESIII analyzed e+e“ 
rjJ/'ip |38f|. Statistically signihcant 7 signals were observed, and the corresponding Born 
cross-sections were measured. In addition, a search for the e+e“ /^Ij process observed 

no signihcant signals, and upper limits at the 90% C.L. on the Born cross-section were set. 


A comparison of the Born cross-sections cr(e’''e —)■ r]J/ip) in this measurement to previous 


results 


13 . 3^ is shown in Fig. [9] (a), indicating very good agreement. The measured Born 
cross-sections were also compared to those of —)■ 7r~^7r~J/'ip obtained from the Belle 

experiment j^, as shown in Fig. [9] (b). Diherent line shapes can be observed in these two 


processes, indicating that the production mechanism of rjj/xp dihers from that of tt+tt J/'ip 
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FIG. 8: Fit to cr(e+e —>■ oJXco) with a resonance (solid curve), or a phase space term (dot-dashed 
curve). Dots with error bars are the dressed cross sections. The uncertainties are statistical only. 



nJs (GeV) 
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FIG. 9: A comparison of the measured Born cross sections of e^e —to those of the previous 

J/Ip from Belle Q]. In these two plots, 


13l . 139(1 (a), and to those of e"''e 


measurements 

the black square dots and the red star dots are the results of rjj/'ip obtained from BESIII. The 
blue dots are results of r]J/ip (a) and vr+vr” J/'i/’ (b) from Belle. The errors are statistical only for 
Belle’s results, and are final combined uncertainties for BESIII’s results. 


in the vicinity of y/s = 4.1-4.6 GeV. This could indicate the existence of a rich spectrum of 
the Y states in this energy region with different coupling strengths to various decay modes. 


VII. MORE DATA FOR XYZ STUDY 


As shown in previous sections, BESIII has achieved a lot in the study of the XYZ states 
and the conventional charmonium states. However, there are more questions to be answered 
with the currently available data and, more importantly, with the data samples that BESIII 
is able to collect in the next few years. 
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A few topics that need to be studied with more data are listed below. 


• In the X sector: 


— Where are the X(3872) and coming from, resonance decays or contin¬ 

uum production? 

— Can the other X states, such as XYZ(3940), X(3915), X(4140) (Y(4140)), and 
Y(4350) be produced in a similar way? 

— Can the charmonium 2P, 3P states and the S-wave spin-singlet states (3S, 4S, 
5S) be observed in radiative transitions? 

• In the Y/'^ sector: 


— Is the Y (4260) structure a single resonance, or does it have a more complicated 
sub-structure? Is Y (4008) a real resonance? 

— What are the other decay modes of y(4360)? 

— What is hidden in the e’''e“ 7r~^7r~hc line shape? 

— Is the Y (4660) observed in e+e“ —7r+7r“-^(2S') the same as the Y (4630) observed 


in e+e —)■ A+A^ 


40 


— What is the correlation between charm production (e’''e open charm final 
states) and charmonium production? 

— Where is the vector charmonium 3D state? 

— Are there charmonium states between '0(4160) and '0(4415), and/or between 
0(4415) and y(4660)? 

— Can the vector charmonium hybrid state be observed dlP*'^ 


• In the Z sector: 


— Are the Zc states produced from resonance decays or from continuum production? 
- Is there a Z^s state decaying into J/0 or D~D*^ + c.c., D*~D^ + c.c.l 
— Are there more Z^. and Z^s states? 

• In the C sector: 

— Can Zlso(2317) be produced and studied? 

— Can the other excited charmed mesons be produced and studied with high-energy 
data? 


BESIII is going to collect about 3 fb“^ of data in the vicinity of 0(4160) to study the Dg 
decay properties in 2015-2016 run. These data can be used to answer some of the questions 
listed above. A few specihc topics are listed here. 

1. In e^e“ —?■ tt+tt” J/0, there is a dip at around 4.17 GeV and a sharp increase at around 
4.23 GeV The data may help to identify where the turning point is located. 
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2. For the e'''e —)> oJXco process (33|, there are few data points close to the threshold to 
support the claim of a narrow structure. 

3. For the e+e“ K^K~J/'ip process, it seems there is a structure at around 


4.2 GeV 42 


4. In the e^e —)■ 'qj/il: process, the line shape differs from that of e^e —)■ tt+tt J/t/', 
with a peak at around 4.2 GeV, but there are no data points on the left side. 

5. In e’''e“ —)■ r]hc, is there a structure close to the threshold? What does the e'''e“ 
rj' J/ijj line shape look like close to the threshold? 


6. Are the e’''e“ IXcj signals from '0(4160) or V(4260), or something else 29 

7. The search for F-wave charmonium states via '0(4160) decays. 

8. The search for 7Xc(2P), 7XYZ(3940), 7^(3915) from 0(4160) decays. 


More data are always better, but as we need to collect data at many points and the 
collection period is limited, we need to dehne lower limits for the number of data points 
and the luminosity needed at each point to ensure meaningful measurements of the physical 
quantities of interest. 

As we require high-precision cross-sections for all of the open charm modes and the cross- 
sections of some hadronic transition modes to £t the resonant parameters of the vector 
states, we require sufficient data points to give the excitation curve of any of the existing 
vector resonances and possible hidden structures in the full energy range. 

Some special issues need to be considered: (1) the thresholds of all the open charm 
and charmonium -|- hadron final states; (2) possible energy regions where large interference 
effects are expected; (3) energy regions where narrow structures are expected; (4) energy 
regions that were not well explored before; and (5) the beam energy spread. 

From previous calculations and measurements, it is known that the energy spread of 
BEPGII at energy ranges of 3.8-5.0 GeV is around 1.4-2.0 MeV. In this case, we would not 
use energy steps finer than three times the energy spread. That is, unless very necessary, we 
would not take data at energy points less than 5 MeV from an existing data point. Thus, we 
would take data with 5 MeV steps in the energy regions where narrow structures or dramatic 
effects are expected, such as the thresholds for open charm or hidden charm final states, 
and the low mass shoulder of V(4260), where interference effects have been reported j^. 
Otherwise, 10 MeV steps will be used. 

The limit on the data size is set according to the precision of the hadronic transition 
modes, which typically have cross-sections of a few to a few tens of picobarns. From the 
analyses in previous sections, we find that an integrated luminosity of 500 pb“^ is needed 
to reach a reasonably high precision for most modes of interest. To search for the XYZ and 
charmonium states via radiative transition, a data sample of at least 500 pb“^ is also needed 
to reach a 5a level observation of a signal if the production rate is 1 x 10“"^ or higher. A 
detailed MG study of the precision that can be reached or the exact luminosity needed for 
each energy point is necessary, as the background level may be very different at different 
resonant peaks. 
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With the above principles in mind, we propose to collect data samples at abont 60 c.m. 
energies from 4.0 GeV to the maximum energy that BEPCII can reach (currently 4.6 GeV) in 
10 MeV steps. Around 500 pb“^ at each energy will be necessary for a comprehensive study 
of the XYZ and charmonium states. As BESIII has already accumulated about 5 fb“^ (see 
Table H), another 25 fb ^ data should be accumulated. This will take about 5 years at 
BEPGII 0. 

VIII. XYZ STUDY IN FUTURE EXPERIMENTS 

Belle-II jdsl will start collecting data in 2018, and will accumulate 50 ab“^ data at the 
T(4S') peak by 2024. These data samples can be used to study the XYZ and charmonium 
states in many different ways [I|, among which ISR can produce events in the same energy 
range covered by BESIII. Figure fTOl shows the effective luminosity at BEPGII energy in the 
Belle-II data samples. We can see that, for 10 ab“^ Belle-II data, we have about 400-500 pb“^ 
data for every 10 MeV in the range 4-5 GeV, comparable to the data sample proposed at 
BESIII in the previous section. Of course, the ISR analyses have a lower efficiency than in 
direct e^e“ collisions because of the extra ISR photons and the boost given to events along 
the beam direction. Even taking these effects into account, the full Belle-II data sample, 
which corresponds to about 2,000-2,300 pb“^ data for every 10 MeV from 4-5 GeV, will 
result in similar statistics for modes like e’''e“ 7r~^n~ J/ijj. Belle-II has the advantage that 

data at different energies will be accumulated at the same time, making the analysis much 
simpler than at BESIII at 60 data points. In addition, Belle-II can produce events above 
4.6 GeV, which is currently the maximum energy of BEPGII. Possible upgrades to increase 
the maximum c.m. energy of BEPGII will obviously expand the physical possibilities of 
BESIII. 



FIG. 10: Effective luminosity at low energy in the Belle and Belle-II T(4S') data samples. 


The HIEPA project 4^1 being discussed at this workshop will improve the aforementioned 
studies in many aspects, particularly with c.m. energies of up to 7 GeV and luminosity 
improvements of a factor of 100. These will allow a hner scan in the full energy region with 
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more integrated luminosity. This will enable a better understanding of all the studies listed 
in this article. 


IX. CONCLUSION 


With the world’s largest data samples at energies of 3.8-4.6 GeV, the BESIII experiment 
made a signihcant contribution to the study of the charmonium and XYZ states. To further 
strengthen such studies, BESIII may collect more data from 4.0-4.6 GeV (or even higher, 
if possible). These data will be complementary to the Belle-II study, with many other 
production mechanisms. The HIEPA project may enable a systematic understanding of the 
nature of the XYZ and charmonium states. 
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